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ABSTRACT 



Aims. We investigate the standing kink modes of a cylindrical model of coronal loops. The density is stratified along the loop axis 
and changes discontinuously at the surface of the cylinder The periods and mode profiles are studies with their deviation from those 
of the unstratified loops. The aim is to extract information on the density scale heights prevailing in the solar corona. 
Methods. The problem is reduced to solving a single second-order partial differential equation for SB,(r, z), the longitudinal compo- 
nent of the Eulerian perturbations of the magnetic field. This equation, in turn, is separated into two second-order ordinary, differential 
equations in r and z that are, however, connected through a dispersion relation between the frequencies and the longitudinal wave num- 
bers. In the thin tube approximation, the eigensolutions are obtained by a perturbation technique, where the perturbation parameter is 
the density stratification parameter. Otherwise the problem is solved numerically. 

Results. 1) On functional dependencies of the dispersion relation the radial wave number is independent of the longitudinal stratifi- 
cation. 2) We verify the earlier computational finding that the first overtone frequencies increase with increasing stratification and the 
observational finding (from analysis of TRACE data) that the ratio of the first to the fundamental overtone frequency decreases with 
increasing stratification. The method we use to arrive at these conclusions, however, is more analytical than computational, and yet our 
numerical results agree with the earlier results. 3) The mode profiles depart from the sinusoidal mode profiles of the unstratified loops. 
This departure and its dependence on the scale height is obtained, and might serve to determine scale heights once high resolution 
data become available. 
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1. Introduction 

Since the earliest identification of the kink oscillations in coro- 
nal loops by Aschwanden et al. (1999a) and Nakariakov et 
al. (1999), a considerable amount of data has been analyzed 
• by Aschwanden et al. (2002), Schrijver et al. (2002), and 
. ^ Wang & Solanki (2004) using the high-resolution observations 
^ of TRACE, SoHO, Yohkoh, etc. See Aschwanden (2003) and 
$H Nakariakov & Verwichte (2005) for an extended review of ob- 
5^ servations of coronal oscillations. 

Verwichte et al. (2004) report periods, phases, damping 
times, and mode profiles for nine coronal loops. As expected, the 
results differ from those based on simplified theoretical models 
assuming cylindrical geometries, constant cross sections, con- 
stant magnetic fields, constant gravitation, isothermal structures, 
constant densities, no initial flows, etc. 

There are numerous attempts to arrive at reasonably realistic 
models where various effects have been from loop geometry to 
structuring and damping have to be addressed. Only after these 
studies can we conclude what may or may not be important in 
the context of solar magneto-seismology, i.e. for solar coronal 
oscillations. Both Smith et al. (1997) and Van Doorsslare et al. 
(2004) have studied the effect of the loop curvature on the oscil- 
lations frequencies. Bennett et al. (1999) and Erdelyi & Fedun 
(2006) studied the twisted magnetic flux tubes in incompress- 
ible media and compare their body, surface, and hybrid modes 
with those of the untwisted cases. Terra-Homem et al. (2003) 
went on to give a detailed discussion of the frequency shifts 
caused by field-aligned background flows. Nasiri (1992) simu- 
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lated a variable cross section by assuming a long narrow-wedge 
geometry. Ruderman (2003) removed the degeneracy inherent in 
loops of circular cross sections by assuming elliptical cross sec- 
tions. Diaz et al. (2001) studied the fast oscillations in the fine 
structure of prominence fibrils. Erdelyi & Carter (2006) then ob- 
tained a full analytical dispersion relation for the propagation of 
MHD waves in structured magnetic flux tubes embedded within 
a straight vertical magnetic environment. Mikhalyaev & Solovev 
(2005) consider the MHD oscillations of double magnetic flux 
tubes in uniform external fields. De Pontieu et al. (2003a, b) an- 
alyze the mechanism of leakage from the photosphere and the 
chromosphere into the transition regions and the corona. Diaz et 
al. (2004) introduce photospheric line-tying boundary conditions 
to emphasize the rate of leakage in damping of the oscillations. 
This is an extension of the infinite homogenous loops of Edwin 
& Roberts (1983). Mendoza - Briceiio et al. (2004) studied the 
effect of the gravitational stratification, and find a 10 - 20% re- 
duction in damping times of oscillations. 

Andries et al. (2005a, b) calculate damping rates of longi- 
tudinally stratified cylindrical loops and conclude that the ratio 
of the frequency of the first overtone to that of the fundamental 
mode is less than 2, the value for the unstratified loops. They use 
this ratio to estimate the amount of the density-scale height in 
the solar atmosphere. Dymova & Ruderman (2005) reduce the 
MHD equations prevailing in a thin and longitudinally stratified 
magnetic fibrils, into a Sturm-Liouville problem for the eigen- 
values and eigenmodes of the fibril. Erdelyi & Verth (2007) use 
the approach of Dymova & Ruderman (2005) to study the devia- 
tions of the mode profiles, i.e. the eigenfunctions, of the stratified 
loops from the sinusoidal profiles. 
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In this paper we study the kink modes of a longitud 
density-stratified loop. We reduce the MHD equations to a 
equation with variable Alfven speed for the z- component 
magnetic field. The dispersion relation, relating the frequi 
and the longitudinal wave numbers, is similar in form t 
for unstratified loops. In the thin-tube approximation, oi 
proach converge to that of Dymova & Ruderman (2005) 
have rescaled the MHD equations from the beginning to a( 
modate thin tubes. 

Equation of motions, boundary conditions, and the rad 
lutions are dealt with in Sects. |2]&[3] The thin tube appro 
tion is treated in Sec |4] Concluding remarks are given in 
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2. Equations of motions 

A coronal loop, with its ends at the photosphere and with 
atively small curvature ( that is, the radius of curvature 
loop much larger than the loop length ), is idealized as a circu- 
lar cylinder. The cylinder is assumed to have no initial material 
flow, to be pervaded by a uniform magnetic field along its axis, 
B - Bz, and to have negligible gas pressure ( zero-yS approxi- 
mation ). The length and radius of the loop are L and a, respec- 
tively, (see Fig.[T]). The coordinate system is the cylindrical one, 
(r, (f), z). The density is assumed to be 



(1) 



(2) 



p{e,z) = pi{e)f{e,z), r<a 
= Peie)fie,z), r> a 

f(e,z) = exp(--sin^), 

TT L 

where p,(e) and Peie) are the interior and exterior densities at the 
footpoints of the loop, e - L/H, where H indicated the scale 
length ( see bellow ). The density variations for inside and out- 
side of the loop are governed by the same function f(e, z). The 
exponential stratification of the density has been pointed out by 
Aschwanden et al. (1999b) on the basis of their EUV studies in 
30 loops from the SoHO/EIT data. The sinusoidal form of the 
exponent is suggested by Andries et al. (2005b). It ensures the 
symmetry with respect to the midpoint of the loop. See, how- 
ever, Erdelyi & Verth (2007) for alternatives to this sinusoidal 
exponent. 

Here, as in most other works, the stratification of the den- 
sity and its exponential scaling is adopted as an empirical fact. 
Evidently its source is not a gravitational one, because, a) the 
gravitational scale height in coronal conditions far exceeds any 
length scale in solar environments and b) the loops are not, gen- 
erally, oriented along the gravitational field of the sun. 

In the following we consider loops of different scale heights 
but of constant total column masses ju, and p^, independent of e. 
Thus, 



Pe 



Jo 



f(e,z)dz 



L[k(-) ■ 
n 



■Lo(-)], 

71 



(3) 



p/(e) Peie) 

where Iq is the modified Bessel function of the first kind and and 
Lo is the modified Struve function (see Gradshteyn & Ryzbik, 
2000). 



The linearized ideal MHD equations for the Eulerian perturba- 
tion in the velocity and the magnetic fields in a zero-/? plasma 
are 



— = —{(V X (5B) X B (V X B) X (5B}, 

of A-np 



(4) 



Fig. 1. A sketch of the equilibrium model of the flux tube. 
Density varies along the cylinder axis and is symmetric about 
the midpoint. The magnetic field is uniform along the z-axis. 



55B 



V X (5v X B). 



(5) 



Note thatEq. (|5]) ensures the Gaussian law, V. 5B = 0. From Eqs. 
(|4]l and Q one can show, after some straightforward calculations 
that 
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d 6B, 
"d'r~B' 



(6) 
(7) 



where va(z, e) - Bj ^Anpiz, e) is the local Alfven speed and is 
different for inside and outside of the loop. The term Svr/iaj ap- 
pearing in Eqs. ^ and (|6]l is actually the Lagrangian displace- 
ment vector in the loop. The remaining components of S\ and 
<5B are given by 



to Id 

6vj, - r6B, —(rdvr), 6v,-0, 

niB- ~ im or 



6Br 



B, d5vr 
ito dz 



5B. 



B, d6v^ 
iu) dz 



(8) 
(9) 



See, e. g., Karami et al (2002) and Safari et al. (2006) for details. 
Let us also emphasize that the present form of Eq. Q utilizes 
the fact that p and, consequently, va depend only on z- There is 
no radial variation except for a step discontinuity at the surface 
of the tube. 

Eliminating 6vr between Eqs. (|7]l and|6]l gives 



I d d d 
r or or Oz 1 



= 0. 



(10) 



This is a wave equation for 6B~ with the variable speed VA{e,z)- 
We solve it by the separation of variables. Let 6B- - R{r)Z{z). 



For later convenience, we add and subtract the term 



v5(f=0) 



6B, 



to and from Eq. ( fTOl ), and follow the usual procedure for the 
separation of variables. We obtain 



I Id^R 1 dR 
R~r~dr^ ^ 'R~dr 

_ 1 d^Z ; 

~ 'zlj''^ 



U) 



0) 
1 



y2(e,z) v2(e = 0)j 



(11) 



where K^{e) is the constant of separation; and in the absence 
of longitudinal stratification, it reduces to the longitudinal wave 
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number. Equation (fTTT i may now be written as 

(12) 



\t2+-1 ^ + fc^/;(r) = 0, 

\dr'^ r dr r'- 

Z(z) + ^Z(z) = 0, k'^ -i— - Ki. 



(13) 



Both Eqs. (fTZb and ( fTlT l are actually a pair for inside, r < a, and 
outside, r > fl, of the tube. They are to be solved simultaneously 
for R, Z, K~, and u. 

2.1. Boundary conditions 

1. The changes in total pressure should be continuous. On ac- 
count of the zero-yS approximation and constancy of B, this 
reduces to the requirement of the continuity of 6Br. Thus 



^interior ^ ^^n^nor ^ f^j. ^jj 



(14) 

(15) 



2. On account of V ■ = 0, should be continuous ?Ar - a. 
This, gives (Karami et al 2002) 



1 dR" 



dr 



1 dR" 



dr 



(16) 



3. The footpoints, z - & L, are expected to be nodes. This 
imposes the conditions 



Z(z = & L) = 0. 



(17) 



Equations ( fT4l i. ( fT6b . and (fTTl i give four boundary conditions 
for the two second-order differential Eqs. ( fT2] l & ( fTSl l. 

3. Solutions of Eqs. Q2) and (QD 

Interior solutions of Eq. ( fT2b that are regular at r = are 7,„(|fe, |r) 
for A:^ > or Imi\ki\r) for A:^^ < 0. Exterior solutions that decay 
with r oo are Kin(\ke\r). They occur for k^ <Q and are evanes- 
cent waves. 

Imposing the boundary conditions of Eqs. ( fT4l ) and ( fTSI l 

gives 



1 y;(|A:,|fl) 1 /:;(|A:,|fl) 

A:; 7m(l^i|a) K,y,(\ke\a) 



0, 



(18) 



where ' indicates a derivative of a function with respect to its ar- 
gument. The same relation holds for surface waves with Jm{\ki\r) 
replaced by Im{\ki\r). For unstratified thin and thick tubes Eq. 
(fTSl l is analyzed by Edwin & Roberts (1983). Here we study 
Eq. (fTSl l for thin stratified loops by perturbational and numerical 
techniques. 



4. Thin tube approximation 

For a/L «c 1 and m > I, the dispersion relation of Eq. ( fTSl l gives 
\ke\. From the definition of k~ in Eq. ( fT2] l. one then obtains 



-[p,(0)+p.(0)]. 



(19) 



This is the kink oscillation frequency in the presence of strat- 
ification. Edwin & Roberts (1978), Karami et al. (2002), Van 
Doorsselaere et al. (2004), and Diaz et al. (2004) all obtained a 
similar result for co in unstratified flux tubes. Here, however, k~ is 



given by Eq. ( fTTT i. It reduces to the longitudinal wave number in 
the absence of stratification. Substituting Eq. ( fT9] l into Eq.(fT3]l 
(interior with kj > and exterior with k^ < 0) and using the 
boundary condition of Eq. ( fTSl l yield 



+ l^F(e, z)Z(z) = 0, Z = Z'""="°'' = Z"'^™^ (20) 
dz -B 

Pi{e)+Pe(e) . 
F(e,z) = f(e,z). 

Equation ( |20] | is an eigenvalue problem weighted by F{e, z). It is 
the same as that of Dymova & Ruderman (2005) derived, how- 
ever, with a different approach. 

From Eq. ( |20l ) one may write down the following integral 
expression for a? 



b2 J \dZldz?-dz 
47r J F{e,z)\Z{z)\W 



(21) 



Some general properties of aP' can be inferred from Eq. (ISTT i. 
In Fig. [3] F{e,z) is plotted versus z for three values of e. It has 
maxima at footpoints and a minimum at the apex. The larger e, 
the higher the maxima and the lower apex become. This reduces 
the integral in the denominator of Eq. (ISTT i and causes an increase 
in a?. There is even the possibility of the integral, and thereby 
w^, becoming infinity. 

In the following section, Eq. ( |20l l is solved for small amount 
of the density scale heights by perturbation, and for arbitrary 
scale height parameters numerically. 

4.1. Perturbation method 

The scale height parameter is chosen as the perturbation param- 
eter, and all variables and equations are expanded in powers of 
e. Thus, 



„((') 



Z(z) = Z(">(z) + eZ<"(z) + 
F(e,z) =po 



1 2 1 . z\ 

1 -I- e — smn— + 

n Lj 

Equation (l20l l splits into zeroth and first-order components 

dz^ B^ 
d^ZP ^ncf^ 
dz^ ^ B2 B2 

.-^^(A_lsin-U<") 



(22) 
(23) 

(24) 



(25) 



(26) 

Solutions of Eq. ( |25] | for Z^°' and aP^ with boundary conditions 
of Eq. (fTTb are 



In 



, --BiAnpo) 



-1/2 



/ = 1,2,3, 



Zf%) 



In 

— sin —z. 



(27) 
(28) 



where I is the longitudinal mode number, and (Jf^ the kink mode 
frequency in the absence of stratification. The right hand side of 
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Eq. ( |26] l is now known. Multiplying it by Z'"'*, integrating over 
z, and reducing it by Eq. ( [25] l gives the first-order corrections 



,(1) 
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7(1) _ V ^ 7(0) ^ _ ' 



(29) 
(30) 



where 



hv 



-7(0)* 1 • ^^-7(0) > 

I Z,, —sin — Z, flz 
Jo L ' 



4 W'(l + cos Ztt cos I'n) 

Z4 + (_l + //2)2_ 2/2(1 + /'2)- 



(31) 



Equation (IJTT i agrees with the result of Andries et al. (2005a) ( 
see 5„ in their Eqs. (3) and (4)). The Sheffield school maintains 
that a knowledge of the deviations of the amplitude profile of 
stratified loops from those of the unstratified one, Z; - zj"' a: 
6Z*'\ in our notation can give information on the density strat- 
ification of the loop (private communication, see also Erdelyi 
& Verth 2007). In Fig. H z|'^ and Z*'* are plotted as functions 
of z- The first, Z' 



<i) 



"1 ^2 

exhibits two maxima at z/L = 
and one minimum at 1/2. It is zero at 0, 1/3, 2/3, 



1/6 & 5/6 
and 1. The 

second, Zj'', shows two maxima at z/L = 1/8 and 5/8, two min- 
ima at 3/8 and 7/8, and is zero at 0, 1/4, 2/4, 3/4, and 1. The 
maxima of zj' '/zf^ and Z^^^ izf are 0.02 and 0.04, respectively; 
see Table [1] column 2 & 4. From the TRACE data, Aschwanden 
et al. (2002) report a displacement amplitude of 100-8800 km 
at the apex of the coronal loops. For a loop of L =100 Mm, 
e - L/H - 2, the calculated percentages, 0.02-0.04, Max(AZi) 
fall in the range 2-176 km. One should be aware of whether the 
accuracies of observed data allows the detection of such minute 
effects. 

We note that w^'* is positive and tends to zero for Z » 1. The 
ratio of the periods of the fundamental and the first overtone is 



Pi 

Pi 



0)2 

a>i 

2 ^ \ = 2 — < 2. 



1 + 6W*'V<' 



(32) 



Equation ( l32b is a useful tool to estimate the density scale height 
of the loops, see also Roberts (2005). 

4.2. Numerical method 

Using a numerical code based on shooting method, Eq. (|20] | is 
solved for eigenvalues and eigenfunctions. For the unstratified 
loops, where p, and pe are constants, the eigenfrequencies and 
the eigenfunctions are those of Eqs. ( |27] ) and ( l28T l. respectively. 
For a range of < e/j: < 25, we have calculated the fundamental 
and the first overtone kink frequencies cj^ and 0J2, respectively, 
and the ratio 0)2/0^1 = Pi /Pi- The results are plotted in Fig. HI 
As anticipated from Eq. ( |2T] ) and the behavior of F(e, z), both 
frequencies show monotonous increase with increasing e. For 
small e, wi has a steeper slope than W2, but both approach each 
other as e increases. The ratio Pi/Pi begins at 2 for unstratified 
loops, e = 0, and decreases to one at large e. From the TRACE 
data, Verwichte et al. (2004) find the ratio 1.64 and 1.81 for two 
of their observed loops. From Fig. ID e corresponding to these 
ratio are 1.93;7r and 1.07?:, respectively. Assuming typical loop 




Fig. 2. Plot of the first-order perturbations in amplitude profiles 
of the fundamental and the first overtone modes. See Eqs. $2M 

Table 1. Maximum amplitude differences between stratified and 
unstratified loops, Ist-order perturbation ( columns 2 <& 4 ), and 
full numerical calculations ( columns 3 & 5 ). 





Max j 
Perturbation 


AZ, \ 

Numerical 


Max| 
Perturbation 


AZt \ 

Numerical 


2 

5 


0.02 
0.02 


0.017 
0.022 


0.04 
0.04 


0.036 
0.04 



lengths, L - 100 - 250 Mm, the density scale height falls in 
the range of H a: 16-4 1 and 30-74 Mm, respectively. These scale 
heights agree with the finding of Andries et al. (2005a, b). 

The longitudinal part of the eigenfield, Z(z), is plotted in 
Figs. |5]for I -1, 2, 3 and e =0, 2, 5. As e increases, a) the 
eigenprofiles depart further from the sinusoidal profiles of the 
unstratified case, b) the antinodes move towards the footpoints, 
and c) the central antinode gets flattened in the case of odd I. 

The differences between the eigenprofiles of the stratified 
and the unstratified cases, AZ/ - Zi{e,z) - Ziie - 0,z), are plot- 
ted in Fig. |6] for the fundamental and the first overtone modes. 
Expectedly, the difference increases with increasing e. The max- 
ima rise and move towards the footpoints as e increases. For 
example, for e - 2 and 5 ( corresponding to L = 200Mm, 
H - 5Q and 20Mm, say ), the first maximum of AZi is located 
at z =36Mm and 35Mm, respectively, in agreement with Erdelyi 
& Verth (2007). 

Table [U shows Max(AZ//eZ/) for / =1, 2 and e = 2, 5. 
Columns 2 and 4 are from the first order perturbation calcula- 
tions. Column 3 and 5 are from the full numerical analysis. The 
proximity of the two different methods of calculations, even at 
scale heights as large as e = 5, is striking. Aschwanden et al. 
(2002) report a displacement amplitude of 100-8800 km at the 
apex of the kink modes. Combined with the fractional deviations 
of Table[T] one may calculate actual physical deviations of 2-176 
km for e - 2 and 6 - 528km for e = 5. This result also agrees 
with Erdelyi & Verth (2007). The question remains as to whether 
the accuracy of the observed data will allow the detection of such 
small effects. 
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e=0 

e=2jt 

e=10Ji 



0.2 



0.4 0.6 
x=z/L 



0.8 



Fig. 3. F(e,z) versus z/L. e = 0.0 (solid line), e = In (dashed 
line), and e = IOtt (dot dashed line) . As e increases, F(e,z) 
recedes towards the end points and is void around where the mid 
point widens. 
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Fig. 4. The frequencies, wi and W2, and the ratio Pi /P2 = (±>2lu)\ 
versus c/tt. Auxiliary parameters are the tube length =100fl, B - 
100 G, and p,/pi = 0.1. 

5. CONCLUSIONS 

We have studied the MHD oscillations of a vertically stratified 
coronal loop 

- Equation ( |20] |. for the longitudinal component of the waves, 
is the same as those of Dymova & Ruderman (2005). 

- The oscillations frequencies, obtained from Eq. ( |20] l, depend 
on the stratification parameter e. This in turn, makes the ra- 
dial wave numbers, k-, & of Eqs. ( fTlT l and (fTST l e- depen- 
dent. 

- In the thin tube approximation, the eigenfrequencies are ob- 
tained by both perturbational and numerical techniques. 

- The effect of stratification is best understood by the behavior 
of F(e, z), highlighted in Fig. [3] F{e,z) is all positive. But 
at large e it becomes insignificant in broad neighborhood of 
z - L/2. This reduces the denominator in Eq. (ISTT i and lets 
oj^ grow. This in turn results in washing-out finite time mea- 
surements of the phenomena under study. 




Fig. 5. a) The fundamental mode (Z = 1), b) the first overtone 
mode (I - 2), and c) the second overtone modes (/ - 3) versus 
z, for unstratified, e = 0, ( corresponding to // = 00, L=100Mm 
, say ) and stratified, e = 2, 5, ( corresponding to // = 50 & 
20Mm and L = lOOMm, say ). 
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Fig. 6. Amplitude differences of the fundamental and first over- 
tone modes from those of the unstratified loops for e = 2 & 5, 
numerical calculations. 



- The ratio of the periods of the fundamental and the first over- 
tone modes (2 for unstratified loops) decreases markedly and 
approaches 1 with an increasing density-scale height param- 
eter ¥or pelpi = 0.1, e/TT = L/ttH = 1.07, & 1.94, the ratio 
f 1/P2 is 1.81 & 1.64. These are in good agreement with the 
observational data of Verwichte et al. (2004), 1 .8 1 ± 0.25 and 
1.64 + 0.23. The latter are deductions from TRACE observa- 
tions, assuming the same density contrast and scale height 
parameter 

- The eigenfunctions of stratified loops deviate from the sinu- 
soidal profiles of the unstratified ones. Relative deviations 
grow with e and are of are close to a few percent, in general 
(see TablellJ. 
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